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Abstract 

Fission-related phenomena of heavy A hypernuclei are discussed with the 
constraint Skyrme-Hartree-Fock+BCS (SHF+BCS) method, in which a sim- 
ilar Skyrme-type interaction is employed also for the interaction between a 
A particle and a nucleon. Assuming that the A particle adiabatically follows 
the fission motion, we discuss the fission barrier height of 23 aU. We find that 
the fission barrier height increases slightly when the A particle occupies the 
lowest level. In this case, the A particle is always attached to the heavier 
fission fragment. This indicates that one may produce heavy neutron- rich 
A hypernuclei through fission, whose weak decay is helpful for the nuclear 
transmutation of long-lived fission products. We also discuss cases where the 
A particle occupies a higher single-particle level. 
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1. Introduction 

Hyperons have attracted much interests since its discovery. The inter- 
action between a hyperon (Y) and a nucleon (JV) is important not only for 
understanding the mechanism of nuclear force but also for the equation of 
state of nuclear matter in astrophysical sites. YN scattering experiments 
have been difficult to perform since the life-time of hyperons in vacuum is so 
short (typically a few 10 -10 seconds). Therefore, the information on the YN 
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interaction has been obtained mainly from the structure of finite hypernuclei, 
which are composed of neutrons, protons and at least one hyperon. There 
have been studied by now many single-A hypernuclei from H to U [lj], and a 
few AA and E hypernuclei have also been produced experimentally 

One of the important questions in hypernuclear physics is to clarify the 
role of hyperon in the structure of hypernuclei. An attractive YN force may 
cause a significant modification in the properties of a finite nucleus when 
a hyperon is added to it. Such impurity effect of hyperon for ^Li, that is, 
the shrinkage of 6 Li core in the hypernucleus, was predicted theoretically 
using the cluster models @, S, S]; which has subsequently been confirmed 
experimentally by Tanida et al. [5]. A recent theoretical work with the 
relativistic mean field (RMF) theory has also shown that the shape of 12 C 
and 28 Si changes from oblate to spherical when a A particle is added to them 



It is intriguing to ask whether the similar effects can be expected for 
heavy fissile nuclei. Interests in the fission of heavy A hypernuclei include i) 
a change of fission barrier height due to the additional A particle and ii) an 
attachment probability of a A particle to each fission fragments after fission 
takes place. The latter has been studied experimentally by Armstrong et al. 
for hypernuclei formed in anti-proton annihilation on 209 Bi and 238 U nuclei 
0]. They reported that the A particle predominantly sticks to the heavier fis- 
sion fragment. Theoretical works on the A attachment probability have been 
also performed in Refs. 0, 0, [h| with a phenomenological Woods-Saxon 
potential and the statistical method, yielding consistent results with the ex- 
periment finding. As for the former interest, to our knowledge, the fission 
barrier height of heavy hypernuclei has been studied neither theoretically nor 
experimentally. 

The aim of this paper is then to clarify the impurity effects of A particle on 
the fission barrier height of heavy A hypernuclei and investigate a possibility 
to produce neutron-rich A hypernuclei, which are difficult to produce directly 
by experimental beams. To this end, we use the (constraint) Skyrme-Hartree- 
Fock+BCS (SHF+BCS) method. This method has been widely used for the 



study of fission I 111 . |12| . Il4 Il5l | and has been applied to hypernuclei as well 



16, 17, 0, 0, 0, 21, 22[. It is also a virtue of this method that one can 
calculate the wave functions of A particle in a self-consistent manner. 

The paper is organized as follows. In section 2, we describe the the- 
oretical framework of the constraint SHF+BCS method. In section 3, we 
illustrate the potential energy surface for fission of heavy hypernuclei and 
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discuss an evolution of A particle motion during fission. In section 4, we 
give the conclusion and discuss possible applications of fission of heavy A 
hypernuclei. 



2. Skyrme-Hartree-Fock method for hypernuclei 

The Skyrme-Hartree-Fock method [13( has been extensively used for sta- 
ble and unstable nuclei. In this approach, the effective NN interaction is de- 
scribed by the Skyrme interaction, which has the form of density-dependent 
zero-range interaction. 



This method has been extended to hypernuclei [16|, in which the A A" and 
ANN interactions are also described by the similar Skyrme-like (^-interaction. 
The two-body AA" interaction is given by, 

v AN {r A ~ r N ) = ^ (1 + x^P a )5(f A - r N ) 

+ ~# (k' 2 5(r A - f N ) + 5{f A - r N )k 2 ) (1) 
+ 4k'5(f A -r N )-k + iW^k'5(f A - r N ) -(axk), 
while the three-body ANN interaction is, 

vann^a,^,^) = t$5(r A -fi)5(f A -f 2 ). (2) 

— * — * 

Here, k and k' are the derivative operators acting on the right and the left 
hand sides, respectively. The parameters in Eqs. ([T|) and ([2]) have been 
adjusted to fit the experimental binding energies of A hypernuclei (l6l . 
With these effective interactions, the total energy E of hypernucleus is given 
by, 

E = E N + E A (3) 
E N = [ H N (r)df+E p / v 

(4) 



E A = J H A (r)dr, 



where E N , E A , and E^ n are the energy of the core nucleus, that of the A 
particle, and the pairing energy, respectively. The energy density for the 
core nucleus, H^ir), is standard and its explicit form can be found e.g., in 



3 



Refs. 13, ll4(- The energy density for the A particle, H\(r), is given by, 

h 2 



H\(r) 



TA + *o (l + \4)PPA + \ (*1 + 4) (TA P + T P a) 



2m a 

+ J (3^ - 4) (Vp • Vp A ) + V A (Vp • J A + V ■ J) (5) 

+ 4*3 Pa (p 2 + 2p„ p p ^ 



where p = p(r) = p p (f) + p„(f), r = r(f) = T p (r) + r n (f), and J = J(r) = 
J p (f) + J n (r) are the number, the kinetic energy, and the spin-current densi- 
ties for the core nucleus, respectively. p\ = p\(r), ta = T A (r), and Ja = J\{f) 
are the same quantities, but for the A particle. 

The single-particle wave functions <p{f) are obtained by minimizing the 
total energy E. This leads to the SHF equations for proton (q = p) and 
neutron (q = n) given by 

("V • + U q {r) + t#(f)) 9 (f) = e q <f> q (^ (6) 

where m*(r) is the nucleon effective mass defined as, 

Uq(r) is the mean-field potential due to the AN interaction given by, 
= 4 (l + 7^0) PtSf) + + 4)r A (f) - 1(3^ - 4)V 2 pA(r) 

- \< (V/ A (r)) + ^<Vp A (f)(-z)(V x a) + ^p A (r) (p(f) + p,(r 

(8) 

Similarly, the wave function for the A particle is obtained with the equation, 

h 2 



■ V *2^ V + l7 A(*0)oA(n =r A o A (r). (9) 



where m\(r) are the A effective mass defined as, 



The mean-field potential U£ for the A particle is given by, 
U»(r) =^(l + ±A p(f) + + 4)r(f) - i(3tf - 4)V 2 p(f) 



+ ^<Vp(0(-z)(V x a) - ^<VJ(f) + ^(p 2 (f) + 2p n (f)p p (? 

(11) 

We define the quadrupole operator as, 



i€q 

and the octupole operator as, 



Ql = E r ' F 3o(^), (13) 



i&q 



respectively (q = p,n,A). Using the self-consistent solution of the SHF 
equations, the multipole moments are calculated as, 



Ql = V^jf / drr 2 Y 20 (6)p q (r) } (14) 
Ql = J drr 3 Y 30 (6)p q (r), (15) 



for the quadrupole and octupole, respectively. 



We perform our calculations by modifying the computer code SKY AX [23, 



24j. This code solves the SHF equations in the coordinate space assuming 
axial symmetry for the nuclear shape. The potential energy surface for fission 
is obtained by constraining on the quadrupole moment, while the octupole 
moment is optimized by adding a small octupole moment to the initial wave 
functions. We take the mesh size of Ar = Az = 1.0 fm with boundary 
conditions of vanishing wave functions at r = 20 fm, z = —25 fm, and z=25 
fm, where r = \J x 2 + y 2 . 

The pairing correlation among nucleons is treated in the BCS approxi- 
mation. We use the density-dependent contact interaction, 

30 = < ( x - 77) - r ^)> ( 16 ) 
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for the pairing interaction. In the code, a smooth cut-off function 



25 



is introduced for the pairing active space. Here, A is the Fermi energy, and 
AE is determined so that 

iVact = J2 /* = N * + l -^ N T ( 18 ) 

k 

with = AE/10, N q being the number of proton or neutron. The pairing 
energy is given by, 

£,pair _ ^2 UkVk A k f k , (19) 

kdp.n 

where u k and v k are the uv factor in the BCS approximation, and A^ is the 
pairing gap. 

3. Results 

Heavy A hypernuclei, which are produced by (K~,tt~), (tt + ,K + ), or 
(e,e'K + ) reactions, are initially at an excited state. Then the fission channel 
competes with several other decay channels. Three fission processes can be 
considered: (I) the A particle itself decays via non-mesonic decay followed 
by fission of the remaining nucleus, (II) fission occurs after the A particle 
de-excites to the lowest single-particle state emitting several 7-rays, (HI) fis- 
sion occurs while the A particle is at an excited single-particle level. Among 
them, the process (I) corresponds to a statistical fission of an excited nucleus 
and we do not consider it in this work. 

We calculate the potential energy curve of 23 ^U as a function of Q 2 = 
Q\ + Q2 + Q21 an d compare it with that of 238 U. We set the pairing strength 
parameters to be Vq = 1410.0 MeV-fm 3 for proton and Vq = 910.0 MeV-fm 3 
for neutron so as to reproduce the empirical pairing gap of 238 U, that is, A n = 
0.674 MeV and A p = 1.168 MeV obtained with the three-point formula for 



nuclear mass. We adopt the parameter set SkM* Jill for the NN interaction, 
while the parameter No. 4 in Table. I of Ref. 17J] for the AA^ interaction, 
which is optimized to 20 ^Pb, 

4 = -253.500 MeV-fm 3 , x£ = -0.216, = 76.920 MeV-fm 5 , 

4 = 28.080 MeV-fm 5 , t£ = .000 MeV-fm 6 , Wq = 0.000 MeV-fm 5 . 
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m tj 


m t j 


B f (inner) (MeV) 


8.20 


8.47 


B f (outer) (MeV) 


6.60 


7.42 



Table 1: The height of the inner and outer fission barriers for the 238 U and j^U nuclei 
when the A particle occupies the lowest single-particle state during fission. 



We have confirmed that the potential energy curve of 23 A U does not alter 



significantly even when other parameter sets for the AN interaction [16|, [17 
are employed. 

3.1. A at the lowest single-particle level 

Let us first discuss the fission of hypernucleus 23 |U at the ground state. 
For this purpose, we assume that a A particle follows adiabatically the fission 
motion. That is, the A particle is at the lowest single-particle level at every 
instant. Fig. [T] shows the fission barrier for 238 U (the dashed line) and 
23 |U (the solid line) as a function of total quadrupole moment Q 2 . Those 
curves are shifted so that the ground state configuration has zero energy. The 
deformation parameter for the ground state does not change significantly by 
adding a A particle. On the other hand, the fission barrier height of 23 A U 
slightly increases as compared to that of 238 U. This suggests that 238 TJ may 
become more stable by adding a A particle, although one needs to evaluate 
also the effect of A particle on the mass inertia to reach a definite conclusion. 
We list the height of the inner and outer fission barriers in Tab. HJ The inner 
barrier height increases by about 0.27 MeV, while the outer barrier by about 
0.82 MeV due to the additional A particle in the 23 |U nucleus. 

The degree of increase of the fission barrier height is primarily determined 
by the energy of A particle, E\. The fission barrier height is defined as the 
energy difference at the saddle point configuration (s.p.) and the ground 
state (g.s.). This reads 

B f = E{s.p.)-E{g.s.) 

= (^(s.p.) - E N (g.a.f) + (£ A (s.p.) - £7 A (g.s.)). (20) 

The first and second terms in Eq. (1201 correspond to the energy difference for 
the core nucleus and the A particle, respectively. We plot these quantities 
as a function of Q 2 in Fig. [2J The top panel shows E N — E N (g.s.) for 
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Figure 1: The fission barrier of 238 U (the dotted line) and 23 j\U (the solid line) nuclei 
obtained with the Skyrme-Hartree-Fock method. The A particle is assumed to occupy 
the lowest single-particle state during fission. The energy curves are shifted so that the 
ground state configuration has zero energy. 



238 U (the solid line) and 23 ^U (the dashed line). Notice that the difference 
between these two curves is small, as shown in the middle panel in Fig. [2j 
The difference is in fact about the order of 1% of the barrier height, and 
thus we conclude that the energy of the core nucleus is insensitive to the 
presence of A particle. In contrast, the energy of the A particle, Ea, varies 
more significantly as a function of Q 2 - The bottom panel in Fig. [2] shows 
E\ — Eji(g.s.) as a function of Q 2 . The energy E\ monotonically increases 
with the quadrupole moment Q 2 with respect to the ground state. The 
energy difference is AE^ = E\(s.p.) — E\(g.s.) = 0.25 MeV at the inner 
barrier position, which almost amounts to the increase of the fission barrier 
height {ABf = 0.27 MeV). It is thus evident that the change in the fission 
barrier height is dominantly resulted from the energy of the A particle. 

Fig. [3] shows the density distributions for the core nucleus 238 TJ (the left 
panels) in 23 ^U and the A particle (the right panels) at the ground state 
(Q2 — 27.15 b), the second minimum (Q 2 = 83.90 b), the outer saddle point 
(Qa = 123.39 b), and Q 2 = 200.00 b. At Q 2 = 200 b, the core nucleus is 
separated asymmetrically into two nuclei, where the fragment on the left- 
hand side is heavier than that on the right-hand side (the initial octupole 
moment at the first stage of iteration determines which fragment is heavier). 
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Figure 2: (The top panel:) The energy of the core nucleus En for the 238 U (the solid line) 
and 23 ^U (the dashed line) nuclei as a function of the total quadrupolc moment Q 2 - The 
A particle is assumed to be at the lowest single-particle state. (The middle panel:) The 
difference between the solid and the dashed curves in the upper left panel. (The bottom 
panel:) The energy of the A particle E\ for 23 aU with respect to that for the ground state 
as a function of Qi- 
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We see that the A particle is localized in the region of z < 0, that is, it is 
stuck to the heavier fission fragment. This is a natural consequence of the 
fact that the binding energy for the A particle is larger in the heavier nucleus. 

Notice that the AN attractive interaction can attract more nucleons in 
the heavier fission fragment. In order to see this effect, we plot in Fig. H] 
the difference between the density of the core nucleus 238 U in 23 ^U and the 
density of 238 U in the absence of the A particle at Q 2 = 200 barn. That is, 

Ap(r)=pi TC (r)-p(r), (21) 

where p^, re (r) * s the density of 238 U in 23 ^U and p(r) is the density of 238 U 
without the A particle. As shown in Fig. [3], the left hand side (z < 0) cor- 
responds to the heavier fragment, to which the A particle is stuck. We can 
see that Ap(r) is positive (red) on the left side, and nucleons are actually 
attracted into the heavier fragment by the presence of A particle. If we inte- 
grate the region of z < 0, about 0.54 proton and 0.86 neutrons are attracted 
due to the A particle. 

The motion of the A particle during fission can also be inferred from the 
single-particle levels for the A particle. We plot them in Fig. [5] as a function 
of Qi- The lowest level (the thick solid line) does not cross with other levels 
in the region which we consider, and the interaction between the lowest level 
and the others is small. This fact validates the adiabatic approximation to a 
large extent. Notice that each single-particle level is smoothly connected to 
those of two isolated fission fragments at large Q2. In particular, the lowest 
level is connected to the lowest level of the heavier fragment because the A 
binding energy is larger in the heavier fragment. As a consequence, the A 
particle is always stuck to the heavier fission fragment when it is at the lowest 
single-particle level. 

3.2. A at a higher single-particle level 

Let us next discuss the fission process when a A particle is at a higher 
single-particle level. This corresponds to the case where the hypernucleus 
goes to fission before the A particle produced experimentally, which is at 
first at an excited state, deexcites to the lowest level. For this purpose, 
we select four single-particle levels with K = 1/2 around the neutron Fermi 
energy, K being the projection of the single-particle angular momentum onto 
the symmetry axis. Those four levels are denoted as level A, B, C, and D 
in Fig. 6 by the thick solid lines. As in the previous subsection, we assume 
that A particle adiabatically moves at the level crossing points. 
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Figure 3: The density distribution for the core nucleus (the left panels) and the A particle 
(the right panels) in 23 ^U at the ground state (Q2 = 27.15 b), the second minimum 
(Q 2 = 83.90 b), the outer saddle point (Q 2 = 123.39 b), and Q 2 = 200 b. The A particle 
is put at the lowest single particle level at all the deformations. 
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Figure 4: The difference between the density of the core nucleus 
density of 238 U in the absence of the A particle at Q2 = 200 barn. 



238 U in 23 ^U and the 




-26 — • — ■ — ■ — • — ■ — • — • — ■ — • — ■ — 1 
20 40 60 80 100 120 140 160 180200 220 

Q2 (barn) 



Figure 5: Single-particle levels for the A particle in 23 ^U as a function of Qi- We assume 
that the A particle stays at the lowest level (the thick solid line) at all the deformations. 
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Figure 6: Same as Fig. 5, but for the cases when the A particle is at four different K=l/2 
levels denoted by the thick solid lines. 



The resultant fission barrier curves for each configuration are shown in 
Fig. [7J The height for the inner and outer barriers is listed in Tab. [2j The 
fission barrier heights for all the cases are lower than that for 238 U shown in 
Fig. [1] (see also Table 1). As we argued in the previous subsection, these 
changes of the barrier height with respect to the barrier height for the core 
nuclues 238 U are intimately connected to the behaviour of A single-particle 
energies. In Fig. [61 one notices that the single-particle energies for the 
occupied levels decrease away from the ground state (Q2 = 27.15 b), leading 
to the decrease of fission barrier height. 

In order to see to which fragment the A particle is stuck when it is at 
an excited level, we plot the A densities for the levels A, B, C, and D at 
Q2 = 200 b in Fig. [HJ We have confirmed that the density distribution for 
the core nucleus is almost the same as that with A at the lowest level shown 
in Fig. 3. In Fig. [BJ one sees that the A particle is attracted by the heavier 
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Figure 7: Fission barrier for the 23 ^U nucleus as a function of the total quadrupole moment 
Q2 when the A particle occupies the four different single-particle levels denoted by the thick 
solid line in Fig. 5. 





(A) 


(B) 


(G) 


(D) 


B f (inner) (MeV) 


7.12 


6.87 


7.31 


8.16 


B f (outer) (MeV) 


5.61 


4.25 


3.95 


3.23 



Table 2: The height of the inner and outer fission barriers for the 23 ^U nucleus for the A 
particle configurations shown in Fig. 6. 
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Figure 8: The density distribution for the A particle at Q2 — 200 b for the configurations 
A, B, C, and D. 
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fragment (z < 0) for the level A, similar to the case of A particle at the 
lowest single-particle level. On the other hand, the A particle moves to the 
lighter fragment (z > 0) in the case of level B. For the levels C and D, the 
density distribution is not localized and we cannot say to which fragment the 
A particle moves. 

Fig. shows the A mean- field potential U^(r) (the solid line) and the A 
wave functions </>a(?", z) (the dashed line) at Q 2 = 200 barn for the lowest, 
A, B, C, and D levels. These are plotted for three different values of r, that 
is, r=0 fm, 2 fm, and 4 fm. The localization of the A wave function for the 
lowest level, the levels A and B are clearly seen. For the levels C and D, 
the wave function is expected to be localized more clearly at larger values 
of quadrupole moment Q2, after a few more level crossings, although we are 
unable to check it numerically because of a difficulty of numerical calculation 
at large deformations. 

4. Conclusion 

We have calculated the fission barrier curve for 23 fU hypernucleus with 
the constraint Skyrme-Hartree-Fock + BCS method and compared it with 
that for 238 U. For this purpose, we have assumed that the A particle adia- 
batically follows the fission process. We found that the fission barrier height 
increases when a A particle is added to the lowest single-particle level of 23 ^U. 
We argued that this is caused because the the lowest A single-particle energy 
increases towards the saddle point. On the other hand, we have confirmed 
that, when the single-particle energy decreases at the saddle point as com- 
pared to the ground state, the fission barrier height with that configuration 
decreases. We also discussed the effect of A particle on the mass partition 
of fission fragments. We have shown that the mass of the fragment to which 
the A particle is attached increases due to the attractive interaction between 
A and nucleons. 

The adiabatic approximation which we employed in this paper would be 
justified when the A particle occupies the lowest single-particle level. For 
a A particle at an excited state, on the other hand, its validity depends 
on how fast the fission takes place. When the fission takes place rapidly, 
the A particle may diabatically follow the fission process. In order to take 
into account the deviation from the adiabatic approximation, a dynamical 
calculation for fission with the Landau-Zener transition at level crossings is 
called for. 
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Figure 9: The mean-field potential for the A particle, Uj^(f), (the solid line) and the A 
wave function, (f>\ (r, z), (the dashed line) at r—0, 2, and 4 fm and at the total quadrupolc 
moment of Q2 = 200 barn for the lowest(the left bottom panel), A, B, C, and D levels. 
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Lastly, we would like to point some possible applications of the fission 
of heavy A hypernuclei. One is a possibility to produce heavy neutron-rich 
A hypernuclei. A light neutron-rich hyp ernucleus \°Li has been successfully 



produced via (tt~,K + ) reaction [26|, |27[. This reaction converts two protons 
to a A and a neutron. It is therefore difficult with this method to produce a 
hypernucleus far from the stability line in the heavy mass region. In contrast, 
since the A particle eventually emerges in one of the fission fragments, the fis- 
sion of heavy A hypernuclei may open up a novel method to produce a heavy 
neutron-rich A hypernucleus. Notice that heavy neutron-rich nuclei (without 
hyperon) has been produced via in-flight fission of 238 U at new generation 
radioactive isotope beam facilities 28|, |29] . Another application of fission of 
heavy A hypernuclei is the nuclear transmutation. If a A particle remains at 
a fission fragment, it eventually decays by weak interaction (predominantly 
by non-mesonic decay). The weak decay of A induces fission, and/or neutron 
and proton emissions of the fission fragment which the A particle originally 
sticks to|30, |3l|. It is an interesting future question whether such processes 



can be utilized to transmute long-lived radioactive wasts produced at nuclear 
power plants. 
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